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Abstract. We present a study of two spectral lines, Fe 1 6173 A and Ni 1 6768 A, that were candidates 
to be used in the Helioseismic and Magnetic Imager (HMI) for observing Doppler velocity and the 
vector magnetic field. The line profiles were studied using the Mt. Wilson Observatory, the Advanced 
Stokes Polarimeter and the Kitt Peak McMath telescope and one meter Fourier transform spectrom- 
eter atlas. Both Fei and Nil profiles have clean continua and no blends that threaten instrument 
performance. The Fei line is 2% deeper, 15% narrower and has a 6% smaller equivalent width than 
the Ni I line. 

The potential of each spectral line to recover pre-assigned solar conditions is tested using a least- 
squares minimization technique to fit Milne-Eddington models to tens of thousands of line profiles 
that have been sampled at five spectral positions across the line. Overall, the Fei line has a better 
performance than the Ni I line for vector magnetic field retrieval. Specifically, the Fe I line is able to 
determine field strength, longitudinal and transverse flux four times more accurately than the Ni I line 
in active regions. Inclination and azimuthal angles can be recovered to ~2° above 600 Mx/cm 2 for 
Fei and above 1000 Mx/cm 2 for Nil. Therefore, the Fei line better determines the magnetic field 
orientation in plage, whereas both lines provide good orientation determination in penumbra and 
umbra. We selected the Fe I spectral line for use in HMI due to its better performance for magnetic 
diagnostics while not sacrificing velocity information. The one exception to the better performance 
of the Fe I line is when high field strengths combine with high velocities to move the spectral line 
beyond effective sampling range. The higher g e ff of Fei means that its useful range of velocity 
values in regions of strong magnetic field is smaller than Ni I. 

Keywords: Instrument: spectral lines, Sun: velocity and magnetic field 



1. Introduction 

The Helioseismic and Magnetic Imager (HMI) needs to accurately measure Doppler 
velocity and vector magnetic field with limited spectral information, sampling 
roughly half a dozen wavelengths across the spectral line. HMI is one of a suite 
of instruments to be included on-board the Solar Dynamics Observatory (SDO) 
with launch currently scheduled for August of 2008. The spacecraft is expected to 
have a large velocity range of ±4000 m/s. 

© 2008 Springer Science + Business Media. Printed in the USA. 
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The first requirement for choosing a spectral line is that it contain a clean contin- 
uum with no blends and no near-by lines. Helioseismology requires that a spectral 
line be narrow and deep since a steep dl/dX ensures greater sensitivity to small 
Doppler shifts. The Ni I 6768 A line satisfied these requirements and was selected 
as the spectral line for the Michelson Doppler Imager (MDI) and Global Oscillation 
Network Group (GONG) instruments. Vector magnetic field measurements benefit 
from a high Lande factor, g e ff, and a simple Zeeman splitting geometry. It is 
also important to minimize the number of blends that become apparent in sunspot 
umbrae where the lower temperatures allow for molecular absorption. 

An understanding of the center-to-limb variations of the selected spectral line is 
important because the accuracy of "look-up" algorithms (such as used by MDI) are 
based on line depth and wing-slope parameters. In addition, it is desirable for a line 
to be relatively insensitive to thermodynamic changes so that heights of formation 
do not drastically change due to moderate perturbations of temperature and density. 
Of course, optimum instrumental transmission at the selected wavelength is also a 
requirement. 

The potential of each spectral line to recover pre-assigned solar conditions is 
tested using a least-squares minimization technique to fit Milne-Eddington models 
to tens of thousands of line profiles that have been sampled at five spectral positions 
across the line. See Figure 1 for an example of filter profiles applied at five positions 
across the Ni I line. Similar research was carried out by Graham et al. , (2002) to 
determine the minimum number of spectral filters necessary to recover the vector 
magnetic field parameters. The work contained in this paper builds on that premise 
insofar as we examine the differences inherent in two spectral lines in their capacity 
to recover the vector magnetic field with filter-polarimetric measurements. 



2. General Information of Spectral Lines 

Five papers discuss the solar spectral Fel 6173 A line (Stenflo and Lindegren, 
1977; Auer et al, 1977; Simmons and Blackwell, 1982; Solanki and Stenflo, 1985; 
Landi degl'Innocenti, 1982). Two papers detail the Nil 6768 A line (Jones, 1989; 
Bruls, 1993). The spectral line central wavelength, effective Lande factor and exci- 
tation potential are noted in these papers and can be found in the first three columns 
of Table 1 . Also found in Table 1 are the line depths and full-width half maxima 
values as seen in the disk center quiet sun Fourier transform spectrometer (FTS) 
atlas. The FTS data were used because the system does not suffer from contamina- 
tion by scattered light. Line depth and width values for the spectral lines may vary 
slightly when observed with different optical systems. 

The heights of line formation are also found in the last two columns of Table 
1. They are estimated using Maltby-M umbral model (Maltby et al, 1986) and 
VAL-C model (Vernazza et al, 1973, 1976, 1981) under non-LTE assumption. 
Stokes profiles were simulated using a non-LTE numerical radiative transfer code 
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(Uitenbroek, 2001) based on the multilevel accelerated lambda iteration (MALI) 
formalism of Rybicki and Hummer (1991). All calculations were made in a one- 
dimensional plane-parallel geometry. The model atoms used were a 52-level atom 
used for Fe calculations and a 25-level atom used for Ni. The calculated results 
are listed in the last two columns of Table I with the first row being the VAL- 
C model height for core and continuum and the second row being the Maltby-M 
model height. 

Table I. Parameters of Fe I 6173 A and Ni I 6768 A 



Wavelength 


geff 


Excit. Pot. 


Depth 


FWHM 


H (core) 


H (cont.) 


(A) 




(eV) 




(A) 


(km) 


(km) 


Fel 6173.34 


2.499 


2.22 


0.66 


0.102 


302 


16 












269 


21 


Nil 6767.68 


1.426 


1.83 


0.64 


0.116 


288 


18 












291 


26 



The calculations of Nil height of formation agree well with other literature 
sources, but the Fel heights differ from Bruls (1993). The Bruls values for Fel 
are h = 238 (core) in the quiet Sun using the VAL-C model and h = 386 (core) in 
umbrae using the Maltby-M model. 

2.1. Line Profiles and Transmission in Sunspot Umbra and in 
Quiet Sun 

Line profiles of Fel 6173 A and Nil 6768 A in sunspot umbra and in quiet Sun 
scanned by the Kitt Peak McMath telescope and one meter Fourier transform spec- 
trometer are available at ftp : //argo . tuc . noao . edu/pub/atlas. Due to 
page limitations, we can not reproduce this data and instead, discuss only the 
important features. 

The profiles of these lines in the quiet sun are relatively clean. Ni I has no nearby 
lines or blends. Fel has a blend 0.6 A away from its line center identified as La II, 
which is at 6172.72 A. This La II blend is the transition of 3F3 - 3F2 where the 
lower level excitation potential is 0.13 eV and the g e jf is 1.5. There are many 
Zeeman components (15-plet). It does not simply split, it broadens as shown in 
the photographic sunspot atlas. Simulations have been carried out to estimate the 
La II blends influence on Doppler velocity measurements. In a worst case scenario, 
results show this blend introduces systematic errors of 1 or 2 m/s if its presence 
is not taken into account. Of course, the HMI look-up algorithm for Doppler mea- 
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surements will account for the existence and behavior of the La II blend in order to 
minimize errors. 




■0.4 -0.2 



0.0 
A 



0.2 0.4 



Figure 1. One possible configuration of the HMI filter set overlayed with the absorption line for Ni I 
6768 A (solid line). The central filter is 90 mA FWHM and the spacing is 76 mA between filter 
centers. 



Since HMI will make vector magnetic field observations, the behavior of these 
lines in umbrae is of interest. It is important to realize that almost all umbral profiles 
are contaminated by molecular blends because the lower temperature of the sunspot 
allows for such spectral absorption not present in the quiet Sun spectra. The Fel 
umbral profile shows an obvious blend in the blue wing. This blend is invisible in 
the quiet Sun spectra. This blend is Eu II at 6173.0 A and there are other two blends 
near this line. Ni I observations in umbra have shown an obvious blend suggested 
to be TiO and two other blends nearby. 

2.2. Trade-off Between Lande Factor and Velocity Range 

There is a trade-off between a spectral line with higher g e fj that performs better 
for magnetic field diagnostics and a line with a lower g e f f for which the velocity 
calculation algorithm performs better. Having a line with a higher g e f t means that 
the Zeeman splitting is greater and the useful velocity range is smaller in the pres- 
ence of a strong magnetic field. When the velocity algorithm no longer responds 
linearly to the Doppler shift due to one wing of the line profile having moved out of 
the spectral sampling range, this is called "saturation". (Note, this is distinctly dif- 
ferent from the saturation experienced using traditional magnetograph algorithms 
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when the Stokes V lobes stop growing but begin separating with increasing field 
strength.) 

A shift in the line profile of 0. 140 A causes one wing of the profile to move out 
of sampling range. (This value is slightly less than twice the filter spacing of 76 
mA. See Figure 1 for the configuration of HMI filters.) We use the Doppler shift 
formula of 5X = \ res t x v/cto determine that Fe I will saturate at approximately 
6.8 km/s while Ni I will saturate at approximately 6.2 km/s due to the Doppler shift 
alone. If the spectral line components are shifted solely due to the Zeeman effect, 
we use the formula 5X = 4.67 x 10~ 13 g e f / A 2 B where 5X and A are in A, g e ff 
is the effective Lande factor and B is in gauss, to determine that a magnetic field 
of approximately 3.1 kGauss will cause Fel to saturate whereas a 4.4 kGauss will 
cause Nil to saturate. Note that Fel has a g e ff of 2.5 while Nil has a g e ff of 
1.4. Of course, when the line is shifted due to both velocity and magnetic field, 
"saturation" will occur at lower values. 



3. Observations of Fe I 6173 A and Observations of Ni 1 6768 A using Mt. 
Wilson Observatory and Advanced Stokes Polarimeter 

The magnetograph at Mt. Wilson Observatory has sufficient resolution and sensi- 
tivity to scan profiles of the spectral lines. The Mt. Wilson system for line profile 
measurements was described in its initial configuration by Ulrich et al. (1991) 
and in its extended configuration by Ulrich et al. (2002). The scanned profile of 
the line Nil in a sunspot umbra shows a pattern in satisfactory agreement with 
that obtained by the Kitt Peak McMath telescope and 1 meter Fourier transform 
spectrometer. The profile also displays all of the fine structure as shown in the Kitt 
Peak's spectrum. 

Fe I line was scanned at Mt. Wilson for quiet sun regions on 2002 December 
4; see the top panel in Figure 2. The Mt. Wilson system uses a moving stage that 
carries a bundle of fiber-optic image reformattors with independent rectangular 
pickups in order to scan the spectrum. The outputs from the pickups are measured 
by ten independent photomultiplier tubes to form independent channels of intensity 
measurement. In 1991 there were two active channels while in 2002/2003 there 
were ten active channels. The spectral range from each pickup is offset from the 
next due to the geometry of the reformattor bundle so that a different range of 
spectrum is scanned by each of the channels. The profiles used here are based on 
the restricted range of spectrum that is covered by all active channels. The scanning 
was alternatively red to blue then blue to red taking a total of 20 minutes. Thus all 
of these profiles are smeared by the five-minute oscillations. The resulting profiles 
are shown in the top panel of Figure 2 for three center-to-limb angles: 0, 45 and 
60 degrees. The line profiles of the Nil were obtained in 1991 at Mt. Wilson for 
the quiet Sun. The resulting profiles are shown in the bottom panel of Figure 2 
for three center-to-limb angles; 0, 45 and 60 degrees. The individual profiles for 
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Figure 3. Fel 6173 A (upper) observed by the ASP. LCP (green) and RCP (red) profiles are shown 
as well as Stokes I (black). The profiles shown here are an average over 100 pixels (pixel size is 
0.6" x 0.37") in the umbra. Nil 6768 A (lower) line in left and right circular polarizations observed 
at Mt. Wilson. 



each sub-scan are not smeared and have a narrower profile. These are a better 
representation of the profiles measured by the HMI filters since each exposure 
is short. These profiles were not used because individually they are noisier than 
the time-averaged profiles and because the exposure duration represented by the 
time-averaged profiles more closely matches the exposure duration for the ASP. 
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Two properties of the Mt. Wilson system have degraded the observed line pro- 
files: 

1. The spectral resolution is about 25 milli- Angstroms due to the finite size of the 
slit and spectrograph. This degradation is corrected by a single pass deconvo- 
lution based on a laser profile in 1991 and based on a sodium cell absorption 
profile in 2002/2003. 

2. The Littrow spectrograph system has internal scattered light which uniformly 
illuminates the entrance apertures of the spectral sampling system. Based on 
a comparison of the very deep Na D lines to observations made with double 
pass spectrograph systems, we find that the scattered light level is about 2% of 
the continuum level. This is an additive offset that should be added to all the 
intensities. 

Stokes parameters I, Q, U, and V were acquired by the Advanced Stokes Po- 
larimeter (ASP) on 2002 March 9 with the Fe I line observing in an active region 
umbra located at S3 W4. The line profiles of Fe I in the umbra in circular polar- 
izations are shown in top panel Figure 3. Also plotted in this Figure is the Stokes 
I in black, as a reference. The profiles are averaged over 100 pixels (pixel size is 
0.6" x 0.37") in this umbra. The left and right-circular polarizations of Nil for a 
sunspot umbra were also obtained at Mt. Wilson, which are shown in lower panel 
of Figure 3. The two blends, the TiO line and CaH line, appeal - to be shifted for the 
left and right circular polarizations, suggesting that the effective Lande factors are 
not equal to zero. 

3.1. Spectral Line Parameters as a Function of Center-to-Limb 
Angle 

The line depth, width, equivalent width, and slopes for the red and blue wings of 
the spectral lines Fe 1 6 173 A and Ni 1 6768 A for the quiet Sun have been measured 
from Mt. Wilson data and can be seen in Figure 4. The line depth and width are 
defined here as the height and width of a Gaussian function which best fits the 
observed line profile. Since measurement for Fe I was only made at three positions, 
0°,45°,60°, in the quiet Sun, we used the Nil data at the same positions for a 
comparison. 

Basically, the line depths for both lines linearly decrease with the cosine of the 
center-to-limb angle, the line widths linearly increase, and the equivalent widths 
slightly increase. The Fel line is generally 2% deeper and 15% narrower than that 
of Ni I. The equivalent width of Fe I is 6% smaller than that of Ni I. The slopes for 
the red wings of the two lines are generally greater than those for the blue wing; 
while the slope of Fe I is greater than that of Ni I. 
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Figure 4. Center-to-limb variation of the lines Fe I 6173 A and Ni I 6768 A in quiet Sun observed at 
Mt. Wilson. The a;-axis is the central meridian distance of the observations as in /i=cos(#). 



3.2. Vector Magnetic Field Retrieval of Nil 6768 A and FeI 6173 
A Using the FeI 6301/6302 A Line Pair as a Reference 

Comparisons of ASP inverted data for the Ni I and Fe I lines show how well the 
lines perform for the purpose of vector magnetometry. ASP inversions are a least 
squares fitting based on the Milne-Eddington solution of the Unno-Rachkovsky 
equations of a plane-parallel magnetized radiative transfer of the Stokes line pro- 
files. Normal ASP operation utilizes the FeI line pair of 6301 A and 6302 A. 
However, another spectral line can be observed simultaneously. On 2002 March 
9, a map of an active region was made by scanning the ASP spectral slit across a 
sunspot near disk center (S3 W4) observing with the FeI 6302 A/6301 A line pair 
and simultaneously scanning the FeI 6173 A line. On 2002 March 10, a map of a 
limb ward active region (S9 E62) was made observing with the FeI 6302 A/6301 
A line pair and Ni I 6768 A line. 
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Figure 5 shows the scatter plots of the values observed with Ni I and Fe I simul- 
taneously with the Fel 6302 and 6301 A line pair. If we assume that the values 
determined from the dual line inversion of Fel 6302 and 6301 Aline pair are the 
"true" values, then it is obvious from Figure 5 that Fel performs better as a vector 
field diagnostic than Nil. This comparison indicates that the higher Lande factor, 
geff, of Fel enables a much better determination of the vector magnetic field and 
filling fraction. 

The increasing uncertainties found in inverted quantities when utilizing a single 
spectral line for inversion, as opposed to two spectral lines, deserves some com- 
ment. When comparing correlation plots of single-line inversion results vs two-line 
results, we see larger discrepancies for single line results with a lower g e ff (1.4 for 
6768) than for a single line with larger g e ff (2.5 for 6173). This highlights the 
necessity for any single line inversion instrument to use a line with a large geff- 



4. Line Performance Using Simulated Profiles 

To estimate line performance, we apply spectral filters (Figure 1) to simulated line 
profiles. The parameters recovered from an ASP observation (discussed in §3) were 
used as input to a ME line profile synthesis code to generate a set of simulated 
Stokes /, Q, U and V profiles for the three spectral lines of Fel 6173 A, Fel 6302 
A and Ni I 6768 A. Filter widths and spacing are scaled by the relative line-center 
wavelengths to obtain an unbiased comparison. We employ only inversion results 
(the ME parameters) for which the same pixel was inverted in the two observed 
wavelength regions providing more than 34000 sets of ME parameters for the 
March 9 data and nearly 22 000 sets of ME parameters for the March 10 data. 

The simulated profiles mimic the observed profiles but lack certain solar affects 
such as molecular blends in the umbra, nearby lines, and asymmetries caused by 
gradients in the solar atmosphere. We have considered that only a fraction, /, of our 
resolution element is occupied by the constant magnetic field (this also accounts 
for scattered light in the instrument). This filling factor, /, the observed contin- 
uum level, and the continuum of the scattered light are all provided by the ASP 
inversion. While this approach mimics the observed profiles in most instances, it 
does not account for other possible distributions of magnetic field inside our reso- 
lution element (Emonet and Cattaneo, 2001; Socas-Navarro and Sanchez Almeida, 
2003). Nor do we explore the effects of atmospheres more complicated than ME. 
However, Westendorp Plaza et al. (1998) shows that ME inversions often provide 
a reasonable mean value of the magnetic field over the line formation region when 
the conditions for the approximation are not satisfied. 

To better test performance over the expected range of HMI observations, we 
assign each pixel a new random velocity in the range ± 6 km/s and random mag- 
netic field orientations. From Figure 5, we see that Nil appears to be observing a 
different magnetic field than Fel. Overall, the Nil parameters have a lower flux 
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Figure 5. Scatter plots of the values observed simultaneously in different lines. The left column 
contains values observed simultaneously with Nil 6768 A and the Fel 6302 A/6301 A line pair on 
2002 March 10. The right column contains values observed simultaneously with Fe I 6173 A and the 
Fel 6302 A/6301 A line pair on 2002 March 9. Field strength (top row), filling fraction (middle row 
- FF)) and inclination (bottom row-Psi) are shown. 
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Figure 6. ASP maps for the continuum at 6302 A. The left map is of NOAA 9856, located at S3 W4, 
and was made at 19:27 UT on 2002 March 9. The right map is of NOAA 9866, located at S9 E62, 
and was made at 18:58 UT on 2002 March 10. The black regions outside of the umbra indicate map 
positions for which inversions were not made. 

density than the Fe I parameters. To remove this bias, we assign the field strength 
and filling factor inverted for Fe I 6302 A/6301 A to each line that we simulate. 
While this gives a more equitable parameter set for comparisons, it may somewhat 
uncorrelate solar connections between magnetic field and temperature (which is 
reflected in the ME thermodynamic parameters). Also, the ASP only inverts when 
the polarization signal is above a threshold of 4 x 10~ 3 of the continuum intensity. 
Therefore, it supplies very few realizations of quiet Sun models. 

Note that we are not creating an inversion technique. Instead, we analyze the 
information content of the measurements -the best possible accuracy that any 
inversion technique would be able to obtain with HMI filter-polarimetric measure- 
ments. To this end, we apply the HMI filters to each ME simulated Stokes profile 
and add anticipated HMI noise based on our assumption of measuring 1.5 x 10 5 
CCD electrons for the quiet Sun continuum. This noise level means that we can 
recover our polarization signal with a sensitivity of 2.2 x 10~ 3 of the intensity 
continuum. This is defined as a simulated observation. 

We employ a weighted Levenberg-Marquardt (LM) least-squares minimization 
to fit Milne-Eddington models to our simulated observation (del Toro Iniesta and 
Ruiz Cobo, 1996). By using the pre-noise input ME parameters as initial guesses, 
we seek to determine how much information was lost by the application of filters 
and noise. The spectral sampling limitation combined with the noise may cause 
the fit to not occur precisely at the input parameters. There may be no well-defined 
minimum. Instead, a range of parameters may all yield the same minimum x 2 
value. The statistical difference between the fitted and the true parameters provides 
a measure of the information content contained in the filtered spectral line. The 
results are verified by setting a new initial guess (original ME parameters + 10 m/s, 
50 G, 2° in each tp and <p, and 0.05 in /) and running a second minimization. The 
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results are essentially unchanged, confirming that our method is not sensitive to 
small changes in initialization. 

We supply the fitting code with both the scattered-light profile and the con- 
tinuum level for the simulated observation which we assume to be previously 
obtained. It is not clear at this time how these quantities will be obtained for HMI. 
It is clear, however, that they are essential for the accurate determination of the 
filling factor and intrinsic magnetic field strength. 



4.1. Line Comparisons 



Before comparing the spectral line performance, note the very different distribu- 
tions of the magnetic field strength, |B|, apparent flux density, F = f x |B| (see 
Figure 7), line depth and width (see Figure 4) for the two active regions. While 
some difference between any two active regions is to be expected, much of the dis- 
tribution difference seen here is due to the different center-to-limb angles at which 
they were observed. Performance comparisons for the simulated observations are 
found in Figure 8. 

In Figure 8, note that all three lines have similar velocity errors for weak fields 
and small input velocities The "saturation" effect discussed in §2.2 is seen in Figure 
8 as the effective velocity range is smaller for the Fe I lines than for the Ni I line 
in the umbra (additionally, the filter set is 10% wider for Nil 6768 A than for Fel 
6173 A). 
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Figure 7. Histograms of ASP inverted ME parameters from spectrally resolved observations of Fe I 
6302 A/6301 A. The top panels are for NOAA 9856 and the lower panels are for NOAA 9866. 
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Figure 8. Spectral-line performance comparison using simulated profiles for Fel 6173 A from 
NOAA 9856 (0), Fel 6302 A from NOAA 9856 (A), Fel 6302 A from NOAA 9866 (box), and 
Nil 6768 A from NOAA 9866 (x). Velocity errors are plotted vs |B| (top left) and velocity (top 
right). Errors of |B| (second row left) and inclination (second row right) are plotted vs apparent flux 
density. Azimuth (third row left) and / (third row right) are plotted vs apparent flux density. Fl vs 
|B| (lower left) and for Ft vs inclination are shown. These results represent the 68% confidence 
level. 
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For magnetic field strength and inclination to the line of sight there is nearly a 
factor of two better performance for the Fe I lines than the Ni I line. Errors decrease 
as the apparent flux density, and therefore polarization signal strength, increases. 
These results reveal an important tradeoff. While Fe I has a better magnetic per- 
formance by a factor of two than Ni I and similar velocity performance for weak 
fields, it will be necessary to "chase the line" during spacecraft motion or suffer a 
less accurate velocity performance in umbrae. 

Note that the Fe I 6302 A performance at disk center is slightly better than for 
the same line at \x = 0.6 and this should be kept in mind when comparing Fel 
6173 A (observed at fj, = 1) and Nil 6768 A (observed at \i = 0.6). Even so, the 
performance difference is very small in most cases and allows for an almost direct 
comparison between all three lines. 

In Figure 8, note that while all lines have azimuth errors less than five degrees 
for strong polarization signals, Nil is noticeably less accurate than Fel for low 
apparent-flux density. / errors decrease for larger polarization signals, as expected. 
For low field strengths, there are relatively low errors for apparent longitudinal flux 
density, Fl = f x |B| x cos tp. In the transition between the weak and strong-field 
regimes, the Stokes V signal becomes more complicated. Once the line is fully 
split, interpretation becomes simpler and, indeed, the errors for Fl decrease for the 
Fe I lines for the strongest fields we examined. 

The errors shown in Figure 8 are the absolute values of the errors. If the absolute 
values were not taken, then the errors should be near zero unless systematic effects 
are present. We do not show a separate plot for the systematic errors, but want 
to mention them. The systematic errors for velocity determination are less than 
1 m/s for all lines in the range of - ±3 km/s. For the extremes of the velocity 
range, the systematic errors underestimate the magnitude of the velocity by about 
3 m/s for all lines.The magnetic field strength as plotted vs flux density (similar to 
second row, left plot of Figure 8 except systematic errors would be plotted about 
the abscissa) shows a systematic overestimation for the Ni I line of 150 G and a 
systematic underestimation for the Fel lines of 50 G at 250 Mx/cm 2 . These errors 
reverse sign above 2 kMx/cm 2 and decrease in magnitude to ±20 G. The filling 
factor errors include a systematic overestimation of / for all lines of about half 
the magnitude of what is plotted in the third row, right column of Figure 8. The 
Ni I line shows no appreciable systematic errors for Fl while the Fe I lines show 
a systematic underestimation by one or two Mx/cm 2 for both very weak and very 
strong fields. These were the only noticeable systematic errors. 

Regarding Ft = f x |B | x sin ip, there is some curious behavior. The Fe I lines, 
especially 6302 A from NOAA 9866, have larger errors at 90° while Ni I errors are 
smaller. Smaller errors are expected as the Stokes Q and U signals will be stronger 
at 90° than at other angles. We hypothesize that this lower performance is due in 
part to the lower signal to noise ratio in general for Stokes Q and U as opposed to 
Stokes V and to the relatively low spectral resolution of HMI. This hypothesis is 
supported by the observation that Fl, which is mostly determined from Stokes V, 
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is better determined by nearly a factor of two. Alternatively, we suggest that Ft 
may be poorly determined for each line for transverse fields of certain magnetic 
field strengths especially near the limb. Since the poor performance is more striking 
for the NOAA 9866 Fel 6302 A data than for the NOAA 9856 Fel 6302 A data, 
this may have something to do with line depths and widths due to the differing 
distributions of magnetic field strengths shown in Figure 7. 

4.2. Magnetograph Mode 

HMI will utilize two CCD cameras. One camera, the "Doppler" camera, will be 
observing only Stokes / and V while running at a higher cadence than the vector 
polarimetry camera. To make performance estimations for this mode of operation, 
we repeat the numerical experiments of the previous subsection only using Stokes 
/ and V information. The LM fits are initialized by the pre-noise input parameters 
altered in the following way: the azimuth is set to zero and not fitted, the inclination 
is set to 2° and |B| is projected onto this new angle. For verification that our results 
reflect the correct x 2 minima, we re-run the fits with the inclination set to a random 
value and |B| is projected correspondingly. This is done because our LM code fits 
separately for |B|, /, and ip instead of simply for Fl. Though not fitting for Fl 
directly may impact the applicability of our technique, the results we obtain from 
these two fittings are essentially unchanged, confirming that we have found the 
correct \ 2 minimum. 

In Figure 9, we find the following differences in line performance for the mag- 
netograph mode compared to when the full Stokes vector is observed. Velocity 
errors are subject to the same "saturation" effect, though now all errors are in- 
creased over their previous levels due to the loss of the information in Stokes Q 
and U. As before, errors for Fl increase until the lines become fully separated. 
It is interesting to note that Fl is determined less precisely without Stokes Q 
and U, even in regions of strong magnetic field. Because of this, Fel lines no 
longer perform better than Ni I for umbral Fl measurements. For very low field 
strengths, the Fel lines perform only about 10 Mx/cm 2 better than the Nil line. 
The performance for Fl is also a stronger function of velocity for the Fel lines 
than for Ni I. The trade-off for magnetograph mode is less pronounced than for the 
vector polarimetry mode. Ni I will better measure velocity in umbra while the Fe I 
will measure Fl better in regions of weak fields. We conclude that the full benefit 
of using the Fe I line only occurs when full vector polarimetry is utilized. 



5. Line Performance Using Observed Profiles 

As in Section 4, we apply simulated spectral filters (Figure 1) to the observed line 
profiles and apply noise such that we can recover polarization with a precision 
of 2.2 x 10~ 3 . The calculations are carried out identically to those of Section 4 
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Figure 9. Performance comparison using only simulated Stokes / and V for Fel 6173 A from 
NOAA 9856 (0), Fe 16302 A from NOAA 9856 (A), Fel 6302 A from NOAA 9866 (box), and 
Nil 6768 A from NOAA 9866 (x). The top panels are velocity errors vs |B| (left), velocity (right). 
The bottom panels are errors of Fl vs |B| (left), velocity (right). The results represent the 68% 
confidence level. 



except that observed profiles are used in place of simulated profiles. Fel 6301 A / 
6302 A data has telluric lines nearby and therefore were not analyzed by applying 
filters. ASP inversion results are used as the initial guess for the LM fit. These 
results using observed profiles, then, provide an indication of the attainable error 
levels for HMI with the proposed filters and polarization precision. 

To estimate errors in the LM fits to the filtered observed profile data, we assume 
that the ASP results are the "true values" and the errors are the difference between 
the LM fits and the ASP inversions. The true performance of the observed data 
lies somewhere between the quoted value and that value minus the uncertainty in 
the ASP inversion. Note that the velocities employed in these calculations are not 
redistributed randomly as done with the simulated profile data. Instead, the veloci- 
ties are the real solar velocities as observed with ASP. In addition to the differences 
between the active regions (see Figure 7), NOAA 9856 for Fe 1 6173 A and NOAA 
9866 for Ni I 6768 A, it must be kept in mind that the Fe I line was observed at disk 
center and the Ni I line was observed at fi = 0.6. 

5.1. FEI6173A 

As seen in Figure 10, the observed profiles performance is poor for apparent flux 
density near 250 Mx/cm 2 . Above 600 Mx/cm 2 , however, field direction (inclination 
and azimuth) can be determined to within a few degrees. For apparent flux density 
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Figure 10. 68% confidence intervals for Fel 6173 A. Errors are shown for the ASP Fel inversion 
data (*), the observed profile data (+) and the simulated profile data (0)- Field strength (upper left), 
inclination (upper right), azimuth (middle left), and filling factor (middle right) are all plotted vs 
apparent flux density. Fl errors are shown vs magnetic field strength (lower left) and Ft errors are 
shown vs inclination (lower right). (Note that the data points for Fl for the observed profiles and 
ASP uncertainty coincide at |B| = 2250 G.) 



above 1500 Mx/cm 2 , magnetic field strength can be determined to within 100 G. 
The observations are less accurate than the simulations, but the performance does 
improve as the polarization signal increases. 

Ft, and Ft have the same dependencies on magnetic field strength and incli- 
nation as discussed for the simulations in Section 4. In general, the accuracy will 
be about 40 Mx/cm 2 for Ft, and 100 Mx/cm 2 for Ft- It appears that the observed 
profile data performance for Fe I should be reasonable for Fl and Ft in general, 
for field directions in all regions with more than 600 Mx/cm 2 , and for magnetic 
field strength in penumbra and umbra. 
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Figure 11. 68% confidence intervals for Nil 6768 A. Errors are shown for the ASP Fel inversion 
data (*), the observed profile data (+) and the simulated profile data (x). Field strength (upper left), 
inclination (upper right), azimuth (middle left), and filling factor (middle right) are all plotted vs 
apparent flux density. Fl errors are shown vs magnetic field strength (lower left) and Ft errors are 
shown vs inclination (lower right). 



As seen in Figure 11, ASP uncertainties for magnetic field strength for this 
line are a factor of four greater for the Fe I line. Oddly, the simulated errors are 
generally as good or better than ASP. This is not the situation for Fel and may 
indicate that uncertainties are coming from the line itself when used for vector 
field determination rather than the spectral filtering of the data. Indeed, the ASP 
ME fits to the observed profiles are of poorer quality for Ni I. There are significant 
asymmetries outside of the umbra and umbral contamination is more troublesome 
than for Fe I. These problems appear to manifest themselves in the magnetic field 
strength observed profiles at 1 kMx/cm 2 and above. 

Nevertheless, field directions are generally recovered as well for this line as for 
Fel, but not until 1 kMx/cm 2 . The difficulties for magnetic field strength are also 
experienced for Ft, and Ft as would be expected. It is interesting to note that the 
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poor Ft performance at 90° exhibits itself for ASP uncertainties for both Ni I and 
Fel. This agrees with our hypothesis from Section 4 that the lower signal-to-noise 
in Stokes Q and U as opposed to Stokes V is partially responsible for these errors. 
From this data, we only predict acceptable performance for Fel 6173 A. 



6. Conclusions 

Both Fel 6173 A and Nil 6768 A profiles have clean continuum and no blends 
which threaten performance (see Figures 2 and 3). The Nil parameters of line 
width, line depth and wing slope show smooth and predictable changes as a func- 
tion of center-to-limb angle (see Figure 4). Ni I has been carefully studied by Bruls 
(1990) who concluded the line was a good choice for helioseismology because 
the line profile is stable and not very sensitive to variation of temperature and 
temperature gradient in the photosphere. Extensive studies of this nature, including 
effects of granulation on the line formation, are yet to be carried out for Fe I. 

Not surprisingly, the higher g e f / Fe I serves as a much better diagnostic tool 
than Ni I to ascertain the inclination angle, field strength and filling factor of the 
magnetic field. This is a robust result as it is obtained repeatedly using a variety 
of analysis techniques. First, the ASP observations of Fe I and Ni I compared to 
the ASP observations of the benchmark Fel 6301 A / 6301 A line pair show Nil 
to have much greater uncertainty in recovering the field strength, filling factor and 
inclination (see Figure 5). 

Second, the simulated line profile data shows Fe I to perform nearly a factor of 
two better than the Ni I line for field strength and inclination determination (see 
Figure 8). Third, results using observed profiles show that Fe I is able to determine 
field strength, longitudinal and transverse flux four times more accurately than Ni I 
in active regions. Inversions using observed profiles show inclination and azimuthal 
errors to be recovered to sa2° degrees above 600 Mx/cm 2 for Fel and above 1000 
Mx/cm 2 for Nil (see Figures 10 and 11). This means, using the Fel line allows 
better magnetic field orientation determination in plage, whereas both lines will 
provide a good orientation determination in penumbra and umbra. From the above- 
mentioned results, we predict acceptable performance for Fel 6173 A for vector 
magnetic field determination. 

In recovering velocities, the performance of the spectral lines from inversion 
methods shows both the Fel and Nil lines perform well for low field strengths 
and low velocities, with errors of roughly 15 m/s. However, the higher g e ff of Fel 
means that its operational range of velocity values in regions of strong magnetic 
field is smaller than Ni I. Fe I performs poorly at high field strengths and high veloc- 
ities where the line has moved beyond the filter sampling range. This implies that in 
order to have accurate velocity measurements in umbrae using the Fe I line, the line 
sampling must be increased by either increasing the number of filter positions or 
chasing the line by changing the filter positions as a function of spacecraft velocity. 
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Simulations using only the Stokes / and V profiles shows that the velocity errors 
increase for both spectral lines in umbra if Stokes Q and U are not considered 
(see Figure 9 as compared to Figure 8) . This result implies that active region 
helioseismology accuracy will improve if the full Stokes vector is used for analysis. 

Using the observed profile performance for Fe I as a good indicator for its per- 
formance on-board the HMI spacecraft, we conclude that it will perform well for 
Fl and Ft in general, for field directions in all regions with more than 600 Mx/cm 2 , 
and for magnetic field strength in penumbra and umbra. 
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